Unconventional magnets in external magnetic fields 

Roderich Moessner 

CN ■ 

Max-Planck-Institut fur the Physik of komplexer Systeme, Dresden, GERMANY 

5h 

Oh E-mail: moessner@pks.mpg.de 

<" 

£f~] • Abstract. This short review surveys phenomena observed when a magnetic field is applied to 

a system of localised spins on a lattice. Its focus is on frustrated magnets in dimension d > 2. 
The interplay of field and entropy is illustrated in the context of their unusual magnetocaloric 
properties, where field-tuned degeneracies assert themselves. Magnetisation plateaux can reveal 
1) ' the physics of fluctuations, with unusual excitations (such as local modes, extended string defects 

»' | , or monopoles) involved in plateau termination. Field-tuning lattice geometry is the final topic, 

where mechanisms for dimensional reduction and conversion between different lattice types are 
discussed. 
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1. Introduction 

The attraction of studying magnets arises from a combination of three factors. Firstly, magnetic 
models often provide the simplest setting for the study of complex cooperative phenomena. 
This has for instance permitted the development of detailed theories of ordering phenomena 
linking microscopies with thermodynamics and even out-of-equilibrium behaviour. Secondly, a 
seemingly unlimited range of experimental compounds [lj in fact realises a wide range of these 
phenomena, the study of which constantly not only throws up new discoveries and but also tests 



o 



^ . existing theoretical understanding. 



Thirdly, a closely connected factor is the availability of a wide variety of probes for studying 
magnets experimentally. Too numerous even to list here exhaustively, it is nonetheless worth 
emphasizing how varied their nature is, ranging from local ones (nuclear magnetic resonance 
or muon spin rotation) to the usual thermodynamic ones. Worthy of particular mention are 
neutrons, which provide information on essentially the full dynamical structure factor in a 
straightforward way, their magnetic coupling being weak enough not to perturb the system 
substantially under usual circumstances. The advent of high resolution X-rays capable of 
providing complementary information about microscopic magnetic structure is one of the current 
exciting experimental developments. 

Magnetic fields provide another very convenient handle for studying and manipulating 
magnets in a controlled and moderately non-invasive way. Often, but by no means always, 
the energy scales are such that readily available laboratory fields, of the order of 10 Tesla, 
are comparable to thermal and exchange energies. Dedicated facilities can provide much higher 
magnetic fields. In particular, pulsed fields routinely achieve strengths well above 50 T. However, 
in this regime, other parameters such as width and shape of the field pulse play an important 
role; for example, for very short pulses, equilibration is not always straightforward. Besides, just 
as impressive as data from successful high-field experiments can be, the aftermath of a failed 
one is also spectacular, but certainly such an experiment no larger qualifies as non-invasive [2]. 



The behaviour of magnets in an applied magnetic field is the subject of this brief review. It 
aims to expose conceptually important phenomena in as simple a setting as possible, using well- 
studied model systems for illustrative purposes. Emphasis is placed on magnets in dimension 
of two or higher, in particular magnets the low-temperature behaviour of which in zero field is 
already sufficiently exotic. This includes in particular classes of frustrated magnets [3] where a 
cooperative paramagnetic or spin liquid regime extends down to very low temperatures [I]. 

We set up the problem in the following section. Section is devoted to magnetisation 
plateaux, including how they are stabilised, and how they terminate. In Section [U we discuss the 
influence a magnetic field can have on the effective geometry of the magnetic system, including 
specific examples of dimensional reduction, and the generation of unusual defects. 



2. The Hamiltonian 

Let us consider spins S, with i = 1...N, residing on a regular lattice. We will consider both 
isotropic (Heisenberg) interactions as well as uniaxial (Ising) ones: 

%u = J z2 &i • Sj 

m 

Hj = Jj2 S i S j , (1) 

m 

where the sum runs over all bonds of the lattice. 

In addition, there is the magnetic field term, which can in its simplest incarnation be written 

as 

N N 

Wfl = = ^ (2) 

where the magnetic moment fi = g^sS related to the spin via the product of g- factor (assumed 
to be isotropic), and the Bohr magneton, \xb- In the following, we will use the reduced field 
h = guBB for notational simplicity. 

In the case of Heisenberg interactions, the field direction defines a privileged axis and 
thus reduces the symmetry of the Hamiltonian from SU(2) to U(l). The situation for Ising 
interactions can be much more complex. In particular, the local z-directions need not be 
identical; as we will discuss in detail below, this can lead to a wide variety of interesting 
phenomena when a magnetic field is applied. 



3. Magnetisation plateaux, and what happens around the edges 

A linear response of the magnetisation, m, to an applied field allows the definition of a 
susceptibility, x = F° r classical magnets at low temperature, one might expect the linear 
regime to be very broad. To see this, let us minimise Hu + T-Lb to find ground states in an 
applied field. By denoting Lij = Si + Sj, and for a lattice of uniform coordination z, we can 
complete the square to write 



n n + n B = -J2 
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From this, it would follow that m 
magnetic field. 

Nonetheless, it is frequently observed that the magnetisation curve as a function of field is 
not at all uniformly linear. The reasons for this can be very varied - the usual suspects are 
thermal or quantum fluctuation, or additional terms in the Hamiltonian omitted in Sec. ([2]). 



3.1. The saturated plateau 

Indeed, the simplest instance of a magnetisation plateau occurs when the field is strong enough 
to saturate the magnetisation: fh = S. Beyond this, it is clearly impossible to increase 
the magnetisation (at any rate without the involvement of some high-energy process such as 
rearranging the crystal-field scheme). Simple though it is, this plateau has some very useful and 
remarkable properties. 

For quantum spins, the saturation plateau is useful as it reveals properties of the Hamiltonian 
in a moderately direct way. This happens because the ground state is a simple, completely 
unentangled state with all spins in the state S z = S: 

\sat} = ®g 1 \S? = S) . (4) 

This is in contrast to the Neel state, for example, which for finite spin S is always dressed with 
quantum fluctuations. Starting from this state, it is then easy to construct excitations. For 
instance, a spin wave with momentum k will have an energy above the field-induced gap given 
by the Fourier transform of the exchange integral, J{k): 

\k) = -^=J2 ^ nS i I sat ) =^ $ I n I $ ~ J &) + const • ( 5 ) 

Here, S~ is the spin lowering operator. 

By measuring the spin wave excitation spectrum in inelastic neutron scattering across the 
Brillouin zone, it is thus possible to extract reliable values of the exchange constants. However, 
this method is still not entirely model-independent as different choices of terms in % can lead to 
similar shapes of J(k), so that one still needs to use educated guesses, e.g. based on quantum 
chemical and exchange pathway considerations, to set up a model Hamiltonian. 

One example of this is provided by a study of CS2CUCI4 [5], a stacked anisotropic triangular 
lattice magnet exhibiting an unusual excitation spectrum in that there appear to be broad 
neutron scattering continua rather than only well-defined magnon peaks (for recent theoretical 
work on this, see [6] and references therein). Besides obtaining values for the leading nearest- 
neighbor exchange, the shape of the dispersion relation together with symmetry considerations 
led to the conclusion that - besides next-nearest and interlayer couplings - there must also be 
a Dzyaloshinskii-Moriya interaction, which arises from spin-orbit coupling. 

In passing we note a different setting where this idea has been developed further. In the case 
of ferromagnets (which of course also have a fluctuation-free fully polarised ground state), one 
can make a fair amount of progress in solving the two-excitation problem, because the single 
spin-flip ones are exactly known. This information is used to infer properties of the ordered state, 
which can be exotic if the bound state of two magnons condenses before the single magnons do 

How does such a plateau terminate as the external field is reduced? The simplest scenario 
for this to occur is via the condensation of the softest excitation. In the case of a square lattice 
antiferromagnet, this would be a magnon at wavevector (tt, 7r), and its condensation leads to the 
symmetry breaking pattern characteristic of a Neel date. As the condensation of this bosonic 
mode breaks the U(l) symmetry of rotations about the axis defined by the field, one can construct 
a formal analogy between this transition and Bose-Einstein condensation. This is explained in 
detail in the review [5J. 

In a frustrated magnet, there is generically no single softest mode. Rather, there can be 
an entire band of minimal energy spinwaves, stretching through the full Brillouin zone. This 
physics is most easily illustrated in two dimensions for the case of the kagome lattice [9], where 
the flat bands are due to local magnon modes around a hexagon. 

Such a local magnon mode is displayed in Fig. ([1]). This mode is local because like in 
an Aharonov-Bohm cage |10j the staggered magnetisation around the hexagon gives rise to a 




Figure 1. Special excitations on 
the kagome lattice. Left: Soft 
fluctuation around a collinear field- 
induced state. Right: Pattern 
of a low-lying local magnon near 
saturation. A magnon on one 
hexagon penalises the presence of 
further magnons on neighbouring 
hexagons. 



cancellation in the exchange field on the triangular sites capping the hexagon. To leading order, 
the number of independent such modes equals the number of unit cells of the lattice: there is 
no way of choosing subsets A, B of the full lattice A made up of the sites contained in disjoint 
sets of hexagons such that A = B and with 4U5 a proper subset of A. Only when A U B = A 
does this become possible, which implies A r mo d es = -Wunitceiis ~~ 1 OH 02] • However, there are 
two further modes winding around the system (for periodic boundary conditions), which gives 
-Anodes = AWitceiis + 1> m keeping with the presence of a flat band and a gapless dispersive band 

At the saturation field, the Zeeman energy cost of reducing S z is exactly offset by the exchange 
energy gain of the staggered transverse magnetic moment. As long as the magnon modes do not 
overlap, their energies are simply additive. By contrast, overlapping modes incur an additional 
energy cost [9], so that there is a degeneracy given by the number of ways of putting non- 
overlapping hexagons onto the kagome lattice. The entropy of the resulting hard-hexagon model 
is known exactly to be S ~ 0.1111 [2j. This entropy can lead to impressive magnetocaloric 
effects, which we will discuss an example of a simpler setting below. Unusual magnetocaloric 
effects due to such anomalous densities of states in frustrated magnets are quite common overall 
- see Reference [15] and citations thereof. 

Let us next consider what happens when the field is lowered further. As the field is turned 
through the saturation field, the energy of the local magnons goes negative, and they thus 
obtain their maximal density. As magnons 'repel' - the magnetisation therefore jumps to 
S — A r max /iV s ites, where N max is the maximal number of non-overlapping hexagons, which equals 
A^sites/9- This magnetisation jump is independent of S, and therfore its relative size vanishes in 
the semiclassical limit, but it is very significant for small S = i! 

3.2. Intermediate plateaux 

The genesis of the saturated plateau is considerably more simple than that of fractional plateaux 
at, say, 1/2 or 1 /3 of the saturated magnetisation. These constitute a subject in their own right, 
and as their origin is very varied, we can only give some illustrative examples here. 

We first discuss fluctuation-induced collinear plateaux. By fluctuation-induced, we mean that, 
whereas energetically a linear m(h) curve would be favourable in a classical magnet, thermal or 
quantum fluctuations instead favour collinear states. However, collinear states are not always 
available given the fixed length of the individual spins. For this reason, the values of fh p for 
which collinear states do exist are stabilised by fluctuations over a finite range of fields around 
hp = rh p /xo- 

3.2.1. Thermal plateaux For thermal fluctuations, this represents a typical instance of classical 
order by disorder [4J. This happens, for instance, in the kagome magnet |16j . the zero-field 
behaviour of which presents a celebrated example of thermal order by disorder in magnetism 



\17\ \18\ IT9| 120], In a nutshell, the frequencies for some fluctuations around collinear states vanish 
in the harmonic approximation (see Fig. [1]). These 'soft modes' in fact typically have a quartic 
energy, E oc 5 4 , so that by equipartition (5 2 ) oc VT, which is much larger that the (5 2 ) oc T for 
quadratic modes at low T. 




Figure 2. Mapping from (down) 
spins on the kagome lattice to (oc- 
cupied) dimers on the honeycomb. 
A dimer rearrangement leaving the 
zero-point energy invariant is also 
shown. The states shown are 
favoured by quantum fluctuations 
in an applied field. 



This means that the fluctuations out of collinear states have a much larger amplitude than 
for non-collinear ones, and hence they lead to a lower free energy for the former. As this free 
energy gain is of entropic origin, it vanishes in the limit of low temperature, and the width of 
the plateau - where the entropic gain outweighs the energetic cost of departing from M oc H - 
vanishes with T. As this plateau corresponds to two spin pointing along the field and one in the 
opposite direction, it occurs at m p = m sat /3. 



3.2.2. Semiclassical plateaux Quantum fluctuations around classical states provide another 
mechanism of stabilising magnetisation plateaux [21] . This follows from the fact that, 
semiclassically, one can use the commutation relationship [S x ,S y ] = ihS z to obtain a 
perturbative treatment of transverse fluctuations. With the classical orientation defining the 
local z-direction, S z = S, the relative transverse fluctuations are parametrically small in 1/5: 



S X Sy 



J • (6) 



Using 1/S as a small parameter yields an expansion of the total energy, schematically 
given by E = E cUssicgX + £ sem iciassicai+ higher order terms, where -Classical is of 0(S 2 ), and 
-^semiclassical ~ £j of O(S), encodes the zero-point energy of the excitations around the 
classical state. 

Unlike the thermal fluctuations discussed above, where the soft modes gave rise to the stability 
of the plateaux, here, all modes contribute. This makes a controlled analytical treatment much 
harder. A long time ago, Henley has argued that the net effect of fluctuations is to induce an 
effective biquadratic exchange, 

«Kff~-^£(&-£i) 2 (7) 

(if) 

favouring collinearity. For the latest analysis of this issue, see [22J. This has been studied 
explicity for the kagome lattice in a field, where a monotonic dependence of zero-point energy 
on collinearity was found for a set of three sublattice states, but of a more complicated form, 
than a simple biquadratic exchange [23] , 



However, this is not the end of the story: there are in fact many different collinear "2-up 
1-down" states with; their number can be determined exactly, because these states map onto a 
soluble problem in classical statistical mechanics: the hexagonal lattice dimer model (Fig. [2]). 
The hexagonal lattice is the lattice defined by the centers of the triangles of the kagome lattice; 
the sites of the latter map onto the bonds of the former. Colouring in the (hexagonal lattice) 
bonds corresponding to a down spin yield a hardcore dimer covering. 

Generalising an approach due to Henley [M] to include magnetic fields, one needs to evaluate 
the zero-point energy given by [25] 



-E'scmicalssical — Tr^j ^— + 3 J 5 a p + hS a p + (S a ^Sj^ — 3S a p) , (8) 

where S a p is the spin connecting shared by triangles a, /3 of the kagome lattice. 

One finds that a uniform "q = 0" state, where all dimers are parallel, minimises -Esemicaissical 
(Fig. [2]). This is somewhat surprising as, up until this result [25J, a "\/3 x \/3" state with an 
enlarged unit cell had always been the fluctuation-favoured state in a wide variety of settings. 
However, related models of quantum fluctuations - such as one for easy-axis spins - have 
yielded a similar ordering pattern since. 

This is not yet the end of the story. There in fact remains a degeneracy - exponential in the 
linear size of the system - involving configurations where an even number of columns of dimers 
have been shifted (Fig. [2j. This results from Henley's gaugelike symmetry [24], and is discussed 
in detail in [23], or, in an analogous zero-field setting [27] . 



3.2.3. Magnetoelastic plateaux A conceptually attractive route to a magnetisation plateau has 
recently been prompted by experiments on the spinel CdCr204 [28} 129] . where the isotropic 
S = 3/2 Cr ions reside on a pyrochlore lattice (Fig. [3]). Here, the plateau is stabilised via 
coupling of the spins to elastic lattice degrees of freedom. In the simplest model, these are 
captured by two additional terms, T-L e encoding the harmonic restoring force for deviations X{ 
of the ions from their equilibrium position, and H me , which models the distance dependence of 
the magnetic exchange via a simple Taylor expansion: 

He = ^ ^ ] ( x i ~ x j) j H-me = J ^ ] (%i ~ x j) S{Sj . (9) 

The phonons can then be straightforwardly integrated out (see e.g. [22]), and one obtains again 
a biquadratic effective exchange, symbolically written as: 

J' 2 ,.^2 

« ~ (SiSj) . (10) 

The further analysis of course does not depend on the origin of this biquadratic exchange. 
In closing, we note that frustrated magnets are inherently susceptible to magnetodistortive 
transitions as their many ground-states include non-uniform patterns of the bond energy Si • Sj , 
and (|10p thus has a nontrivial action in the space of degenerate states to first order in its 
prefactor [30] [31] . 



3.2.4- Commensur -ability plateaux There is a long-standing strand of activity centered around 
the "Shastry-Sutherland" [32] compound SrCu2(B03)2, where various intermediate plateaux 
have been, and are being, discovered experimentally. Their theoretical study is being pursued 
under the heading of supersolid magnetic phases. This field is too broad, and too much in flux, 
to be given justice here, and we refer the interested reader to an upcoming review [33j . 



The basic idea again is to treat excitations as weakly interacting Bosons [8] . At commensurate 
densities p = p/q (where p, q are small integers), a Mott transition for these Bosons can occur, 
the incompressibility of which manifests itself as a magnetisation plateau. 

Incompressible states are of course not exclusive to Bosons - indeed, Fermions in an external 
field exhibit the fractional quantum Hall effect, itself one of the most prominent incompressible 
states. Based on a loose analogy to this physics, Misguich et al. [33] have developed a Fermionic 
theory of the magnetisation plateau. The status of the connection between Fermionic and 
Bosonic theories is at present unclear. 

This concludes our survey of Heisenberg magnetisation plateaux. The next section, which 
deals with the interplay of magnetic field and lattice geometry, we will encounter another plateau, 
which goes along with a reduction of the effective lattice dimensionality. 

4. Tuning lattice geometry via a magnetic field 

To demonstrate the degree of control magnetic fields provide, this section discusses how they can 
be used even to modify the effective lattice geometry. In the following, we show how the lattice 
type can be tuned (pyrochlore to face-centered cubic), how its dimensionality can be modified 
(dimensional reduction from d = 3 to d = 2 or d = 1), and how the nature of the low-energy 
degrees of freedom changes along the way; we will e.g. encounter extended string defects and 
spins fractionalising into monomers/monopoles. 

The mechanism whereby this is achieved involves coupling the field selectively to certain 
subsets of spins. The central ingredient to our example is a classical Ising model in which the 
Ising axes are non-collinear |35| . This can be realised for four spins arranged on the corners of 
a tetrahedron, which has four such non-collinear axes (see inset of Fig. [5|). These are the cubic 
[111] axes, d K (k = 1...4), which all have a mutual angle, 9, with cos# = dk ■ d' k = —1/3. 

Figure 3. The pyrochlore lattice. 
One sublattice (green dots) forms a 
face-centred cubic lattice. Orthog- 
onal chains are shown as red lines 
(solid and dashed). One kagome 
star in a {111} plane is shown 
(dashed blue) , along with a triangle 
belonging to the alternating trian- 
gular plane (dashed green). 

This non-collinearity has as its first striking consequence that a ferromagnet is frustrated, and 
an antiferromagnet is not. Here, by (anti) ferromagnet, one means a positive (negative) Curie 
temperature, O oc zJ, obtained from the inverse susceptibility \ oc T — O at high temperature 
T S> |0|. Introducing a variable Uj = ±1 which encodes if a spin points into or out of the 
tetrahedron along its Ising axis, one finds 

(ij) (ij) (ij) 

Hence, the low-temperature is described by an Ising model with J e g = —J/3, for a 
ferromagnetic J < 0, one thus obtains an antiferromagnetic J e g. For the single tetrahedron, 

— ^ a i a j = ~^ \ J2i=t a ij + it ^ere are thus (tj = 6 ground states, those with two spins 
pointing in, and two out. 




Now, consider applying a magnetic field along one of these axes, so that its projection onto 
the spin on the corresponding 'sublattice' (i = 1, say) is three times that onto the others. A 
magnetic field which is much smaller than J will lift the degeneracy of the six ground states. 
Having fixed the direction of one spin with the field (a = 1, say), exactly one of the three 
remaining spins must have the same value a = 1, and the ground state degeneracy is thus 
halved. 



Next, let us combe these tetrahedra to form a lattice. A particularly productive way of doing 
this is to arrange the tetrahedra so that they share corners (Fig. [3]), with their centers forming 
a diamond lattice. The resulting pyrochlore lattice decomposes into the pinned sublattice, 
which lie on a face-centred cubic lattice. This can be thought of as stacks of triangular layers, 
alternating with kagome layers hosting the spins on the other three sublattices. If interactions 
exist only between nearest neighbours, these layers are decoupled thanks to the inert triangular 
layers-dimensional reduction from d = 3 pyrochlore to d = 2 kagome has occurred. This field- 
tuned dimensionality is also evident in the magnetic correlations |36t [37] 

With each triangle harbouring exactly one spin with a = +1, we retrieve the mapping to a 
dimer model on the hexagonal lattice described above in the context of semiclassical plateaux. 
Hence, the plateau exhibits a finite entropy per spin of the pyrochlore lattice of S ~ 0.08£;£. 

This entropy is different from that of the zero-field case, which has an entropy per spin of 
S p ^log|. This "Pauling entropy" was derived by Linus Pauling in the context of ice, and 
indeed the Ising pyrochlore antiferromagnet is known as spin ice: tetrahedra with two spins 
pointing in and two pointing out are said to obey the ice rule. This model was first discussed 
by Anderson in 1956 [38J, with the corresponding compounds discovered in 1997 [39|, and the 
Pauling entropy measured shortly thereafter |40| . There has been a sustained research activity 
ever since, which we will have to review on a separate occasion. 

How is three-dimensionality re-established when the field strength is lowered at low T? This 
happens via string defects, which are themselves one dimensional. Their extended nature follows 
from a simple consideration: the tetrahedra pointing upwards (and downwards) from a given 
triangular layer comprise all sites of the adjacent kagome layers. As J2i=i °~i = f° r each 
tetrahedron it follows that the magnetisations of kagome and triangular layers are equal and 
opposite; iterating this, one finds that all triangular layers have the same magnetisation. Hence, 
under the constraint J2i=i °~i = everywhere, one can only flip strings containing at least one 
spin in each layer. Such strings need not be straight: they can meander in the kagome layers, 
and the entropy of this meandering in fact sets their density at |36j : 



entropy 
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Figure 4. Schematic behaviour of 
magnetisation and entropy for spin 
ice in a [111] field. The central 
"kagome ice" plateau is terminated 
by string and monomer defects at 
low and high fields, respectively. 
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n str ~ exp [32hJ/9kBT] 



(12) 



The termination of the plateau at high fields occurs by a different mechanism. Now, the 
Zeeman energy h ■ d^uscTi wins over the exchange, and the ground state becomes the unique 
saturated state with u\ = 1 on the triangular layers and cr, = — f on the kagome ones. A 
triangle with all spins a = — 1 , in dimer language, does not have a dimer attached to it. For this 
reason, such a configuration is called a monomer defect - this is now a point-like object. The 
energy difference between this and a tetrahedron obeying the ice rule is given by 

AE = 2\J\/3-h/3 . (13) 

At h = /i sat = 2J, the (ice rule) state is degenerate with the (3-in, 1-out) state and one can 
choose which tetrahedron is in what state. Just as before, this local choice gives rise to a peak 
in the entropy, which here is gigantic: S ~ 0.291/cb, in excess even of the Pauling entropy. 

Neglecting other (higher energy) configurations, it follows from the form of Eq. (|13|) . that 
the partition function, and hence all thermodynamic quantities, near saturation depend only on 
the combination {h — h sa t)/T. Under adiabatic conditions, when S is constant, {h — /i sat )/T will 
also have a fixed value. This is completely analogous to the case of a paramagnet, where Z is a 
function of h/T only. 

Cooling by adiabatic demagnetisation uses this fact by reducing h at fixed S, so that T is 
reduced - in an ideal setting - in lockstep. The same effect occurs here, the difference being that 
h is reduced to /i sat / to effect cooling, and hence cooling in a finite final field becomes possible. 
Indeed, the same argument applies for cooling by adiabatic magnetisation. As h approaches h sa t 
from below, a drop in T also occurs! 

This theory has been put to an experimental test which it has passed only partially, 
namely in an intermediate temperature regime. At low T, the crossover captured Eq. (|13j) was 
found to be replaced by a first-order transition, which terminates in a critical endpoint just like 
in a liquid-gas phase diagram (Fig. [5]). 

This has turned out to signal a set of interesting properties of the monomer excitations. The 
fundamental observation is that the Hamiltonian in Eq. ^) is not complete - indeed, the leading 
interactions between spins are long-ranged dipolar ones [10] • Why these have otherwise little 
influence on the behaviour of spin ice is discussed in Refs. \42\ l4"3l 144] . For our purposes, it 
is sufficient to note that two tetrahedron with (3-in, 1-out) configurations acquire a Coulomb 
interaction 

/xo(^) 2 
4i r 

where /2 is the magnetic moment of the ions, and the distance between the centers of two 
neighbouring tetrahedra. Note that the prefactor contains no - this is a magnetic Coulomb 
interaction, and the monomer excitations interact like magnetic charges - or magnetic monopoles 



Whereas the defect-free kagome ice configuration is effectively two-dimensional, the 
monopoles also interact with each other between the planes. This is thus a system where 
the effective dimensionality in the presence of defects changes. The idea that low dimensionality 
obtains only at low T is of course commonplace. Another famous example in frustrated 
magnetism is provided by Lee et al. [45] . 

It is even possible to generate a dimensional crossover (one-dimensional) chains by applying 
a judiciously chosen magnetic field, namely one in a [110] direction, which has zero projection on 
two sublattices and a finite and equal projection onto the other two; these latter ones get pinned, 
and the former point head-to-tail along chains running perpendicular to the field direction (see 
Fig. E|). The overall direction the spins of a given chain point in is a priori undetermined, and 
one can thus define an effective Ising variable /ip = ±1 which encodes either option for each chain 
r. A detailed numerical study for the dipolar spin ice Hamiltonian has found that neighbouring 
chains order antiferromagnetically, thus establishing three-dimensional order at low T [46j . 
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Figure 5. Liquid- gas phase dia- 
gram of spin ice with long-ranged 
dipolar interactions in a mag- 
netic field. For nearest neigh- 
bour interactions only, the liquid- 
gas phase transition gets replaced 
by a crossover. Inset: non-collinear 
Ising axes (dashed lines) of spins on 
a tetrahedron. 



We conclude this parade of field-induced effective lattice geometries with the face-centered 
cubic lattice, which can be obtained if the field is rotated away from [110] so that three sublattices 
are pinned by it, with the fourth sublattice experiencing a field of a strength to cancel exactly 
the exchange field due to the other three sublattices. Indeed, this is the situation one obtains 
also when one tilts the field at h sa t away from [111] so that one kagome sublattice is pinned more 
weakly than the other two. The undetermined sublattice now forms a face-centred cubic lattice, 
and the spins on it elect to order spontaneously in the presence of long-range interactions |46| . 

5. Conclusion 

We hope this review provides a useful idea of the versatility that magnetic fields offer in the 
study of magnetism. As always, much more has had to be omitted than was included. Obvious 
topics which merit a more detailed treatment include magnetoelastic effects, the role of defects, 
properties of more complex models systems, exotic field-tuned transitions and of course an in- 
depth treatment of quantum dynamics. Regarding the latter, we would like to mention the 
possibility of generating quasiparticles with non-Abelian statistics by applying a magnetic field 
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